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ABSTRACT 


This  paper  Investigates  in  what  manner  the  spectrum  of 
body  waves  radiating  from  point  sources  in  a  multilayered 
medium  ovep  a  homogeneous  half-space  is  different  from  the 
spectrum  of  body  waves  from  the  same  source  in  an  infinite 
medium.  The  effect  of  the  system  consisting  of  a  point 
source  in  a  layered  crust  on  the  spectrum  of  p-waves  ob¬ 
served  at  large  distances  in  the  half  space  is  studied. 
Analytical  expressions  for  the  transfer  function  of  this 
system  are  derived  for  three  types  of  point  sources:  an 
explosive  source,  a  single  couple,  and  a  double  couple  of 
arbitrary  orientation  within  the  crust. 

Preliminary  numerical  computations  for  the  explosive 
source  at  various  depths  in  a  realistic  model  of  the  earth's 

t 

crust  study  the  effect  of;  a)  the  angle  of  incidence  into 
the  homogeneous  half-space,  b)  the  source  depth,  c)  minor 
variations  of  the  crustal  model.  In  the  case  of  an  explo¬ 
sive  source  the  most  influential  parameter  of  the  transfer 
function  is  the  source  depth.  In  shallow  explosions  the  low 
frequency  part  of  the  spectrum  of  body  waves  is  compara¬ 
tively  rejected. 
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Point  Sources  and  the  crust  of  th  earth. 

An  explosive  source,  a  single,  or  4  double 
couple  radiate  body  waves  Into  the  mantle 
with  angle  pf  Incidence  ]f  p. 

Cartesian  (x,y,z)  and  cylindrical  polar 
(T1O4Z  )coordlnate  system  with  numbering 
of  layers,  depth  of  interfaces  and  of 
source . 

Contour  Integration  in  the  complex  -plane. 

Geometry  of  the  single  couple.  The  point 
forces4JL^and-|||^  act  on  the  line  element 

The  Influence  of  the  angle  of  Incidence 
on  the  anqplitude  of  the  transfer  function 
of  the  system  consisting  of  an  explosive 
point  source  in  the  crustal  model  NOCR 
(right  part  cf  figure)  -  at  depth  30  km. 

t  =  75°,  60°,  45°,  30°,  15°  from  top  to 
bottom.  Corresponding  phase  velocities 
are  indicated. 


Anqplltude  of  the  transfer  function  for  the 
explosive  source  at  depths  0.5,  8.0,  30.O  km 
in  the  crustal  model  NOCR.  Angle  of  incidence 
I  =  45°.  Frequency  range  0-0.3  cps. 

Anqplitude  i  T  |  and  phase  4  of  the  transfer 
function  for  the  explosive  so^ce  at  depth 
0.5  tan  in  the  crustal  model  MOCR.  Angle  of 
incidence  J  =  45°.  Frequence  range 
0-4.0  cps. 

Amplitude  m  and  phase  f  of  the  transfer 
function  for  the  explosive  source  at  depth 
8.0  km  in  the  crustal  model  i^OCR,  (See 
Figure  7) 

Congjarlson  of  amplitude  of  the  transfer 
functions  for  explosive  source  at  depth 
8  km,  f  =  30®  in  crustal  model  NOCR,  and 
modified  models  LVCl,  LVC2  (see  tables  I-3). 


The  Transfer  Function  for  P-Waves  for  a  System 
Consisting  of  a  Point  Source  in  a  Layered  Medium 


1.  Introduction. 

Focal  mechanism  studies  using  body  and  surface  wave 
observations  endeavor  to  extract  useful  information  about  the  | 

source  of  an  earthquake  from  the  set  of  signals  observed  at  | 

seismic  stations  distributed  over  the  world.  Not  all  the  * 

information  about  the  source  originally  contained  in  the  sig-  j 

nal  in  the  neighborhood  of  the  source  will  reach  the  stations  j 

at  large  distances  from  the  sox^rce.  There  is,  of  course,  I 

loss  of  information  due  to  geometrical  spreading  and  to  ab-  ^ 

sorption.  But  equally  as  Important,  each  homogeneity  encoun-  | 

'  I 

tered  by  the  signal  along  its  path  from  the  source  to  the  j 

receiver  will  more  or  less  reshape  the  signal.  A  very  strong  j 

deformation  of  the  seismic  signal  occurs  on  its  passage  | 

£ 

through  the  crust  of  the  earth,  since  the  dimensions  of  the  I 

I 

crust  are  of  the  order  of  the  dominant  wave  length  of  the  I 

signal.  j 

For  the  case  of  deep  focus  earthquakes  (source  | 

i 

depth  >  300  km)  only  the  crust  at  the  receiver  has  to  be  j 

accounted  for.  But  as  soon  as  the  source  approaches  the  .  | 

crust.  Interference  between  direct  waves  and  reflected  waves  I 

becomes  a  factor.  In  this  circumstance  the  structure  of  the  I 

earth's  crust  at  the  source  will  also  strongly  deform  the  | 

I 

signal  originally  radiated  from  the  point  source.  Espc-  | 

dally  is  this  the  case  if  the  source  is  located  within  the  | 

K 

I 
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crust.  It  is  the  purpose  of  this  study  ^ '>  Investigate  how 
the  spectrum  of  body  waves  radiating  from  a  point  source  is 
affected  if  this  point  soux’ce  is  placed  within  the  crust  in¬ 
stead  of  in  an  infinite  homogeneous  and  isotropic  medium. 

This  paper  concentrates  on  the  spectrum  of  dllatational  waves 
radiated  from  three  commonly  used  point  source  models  (e^.plo- 
slve  source,  single  couple,  and  double  couple)  into  the 
mantle . 

The  dllatational  displacement  potential  in  the  mantle 
of  the  earth  at  large  distances  fr'>,.i  the  source  will  be  cal¬ 
culate'  .  Numerical  examples  for  the  case  of  the  explosive 
source  in  the  crust  will  be  presented.  For  the  sake  of  clar¬ 
ity  the  assumptions  made  throughout  this  paper  concerning  the 
nature  of  th?  model  consisting  of  a  point  source  in  the  crust 
will  be  listed  here. 

1 . 1  Assumptions  on  the  model  "point  source  in  the  crust.” 

The  situation  is  depicted  in  figure  1.  The  assump¬ 
tions  are; 

1.  We  are  dealing  with  a  point  source,  l.e.  the  dimen¬ 
sions  of  the  source  region  should  be  small  compared 
to  the  distance  to  the  next  boundary.  With  this 
limitation  in  mind  the  point  source  may  be  placed 
at  any  point  in  a  layered  medium. 

2.  The  crust  of  the  earth  is  approximated  by  a  system 
of  (n-l)  homogeneous,  isotropic  layers  with  plane 
parallel  Interfaces.  The  P-,  the  S-wave  velocity, 
the  density,  and  the  layer  thickness  for  the  1 


Source 


J»  V) 

-  a* 

tk  o  , 


V 


-r: 

S  k 


2  *5 

12 


c  a 


S  0: 


V) 

Uj 

o 

Q: 

o 

</) 


»  «**  i 


c 

o 

«  o 


>  ^ 

^  :d 

o  Q: 

0.  o 


3 


layer  will  be  dor.oted  by  C^j^,  (^i>  ^1  ^1 

spectively.  The  top  layer  is  specified  to  have  a 
free  surface.  The  angle  is  the  angle  of  inci¬ 
dence  for  P-waves  in  the  mantle. 

4*  w 

3.  The  mantle,  the  n  layer,  will  be  regarded  to  be 
"relatively  homogeneous."  By  this  we  mean,  quite 
generally  speaking,  that  the  elastic  properties  of 
the  mantle  vary  much  more  gradually  than  they  do 

^  Vi 

in  the  crust.  Thus  we  choose  the  n  layer  to  be 
a  homogeneous  half-space. 

Prom  this  point  on,  the  word  "crust"  will  mean  the  system 
of  n-1  layers  described  under  1.  The  term  "mantle"  or 
"subcrustal  material"  refers  to  the  homogeneous  half-space 
of  our  model. 

Figure  1  .  indicates  furthermore,  how  the  signal 
might  be  followed  through  the  mantle  to  the  base  of  the 
crust  at  the  receiver  station  and  through  that  crust  to  the 
receiver.  In  this  paper  vie  will  not  consider  the  Influence 
of  the  mantle  path  or  of  the  receiver  crust  on  the  signal. 
These  latter  effects  have  been  examined  by  a  number  of 
authors.  Ben-Menahem,  et  al.  (1965)  descrlb'^d  a  method 
for  the  recovery  of  source  parameters  if  the  source  is  lo¬ 
cated  Ir  a  realistic  mantle  model  sufficiently  removed 
from  the  base  to  the  crust.  The  influence  of  the  receiver- 
crust  system  has  been  treated  recently  by  Thinney  (1964), 
Kannon  (1964a,b),  Ben-Menahem,  et  al.  (1965)  using  the 
matrix  formulation  of  Thomson  (1950)  as  developed  by 
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j  Haskell  (1953). 

I  1.2  Previous  Investigations. 

I  Elastic  waves  radiating  from  point  sources  In  a  mul¬ 

ti-layered  medium  overlying  a  homogeneous  half-space  have 
been  Investigated  recently  by  a  num!  er  of  authors.  Hark- 
rlder  (1964)  formulated  the  solution  for  an  explosive  source, 
a  vertical  and  a  horizontal  point  force  In  a  layered  medium 
by  the  Thomson-Hasi^ell  matrix  method  In  the  frequency  domain. 
He  concentrated  on  the  effect  of  the  source  crust  on  surface 
waves  of  Flaylelgh  and  Love  type. 

Bortfeld  (1964)  reported  on  a  solutJ  In  the  time 
domain  by  a  method  of  numerical  interpretation. 

Van  Nostrand  (1964)  presented  computations  of  synthet¬ 
ic  seismograms  at  large  distances  from  point  sources  in  the 
crust.  He  also  works  completely  in  the  time  domain. 

Dunkln  (1965)  applied  integral  transformation  tech¬ 
niques  for  a  refinement  of  the  matrix  formulation.  His 
matrix  formulation  is  especially  useful  for  numerical  cal¬ 
culations  at  high  frequencies.  He  provides  a  review  of  the 
current  Thomson-Haskell  matrix  method  which  gives  a  good 
insight  into  the  purpose  and  significance  of  matrix  methods 
in  problems  of  wave  propagation. 

In  this  paper  we  shall  now  proceed  in  the  following 

steps; 

Section  2:  The  displacement  potentials  for  P-waves  die  to 


a  vertical  and  a  horizontal  point  source  In  the 
crust  will  be  derived  as  integral  representation 


(sections  2.1  and  2.2).  A  large  distance  approx¬ 
imation  will  be  given  in  2.3. 

Section  3:  The  results  for  the  vertical  and  horizontal  point 
force  are  used  to  construct  the  potentials  due  to 
three  types  of  point  sources  commonly  used  in 
focal  mechanism  studies:  an  explosive  source, 
a  single  couple,  and  a  double  couple.  Transfer 
functions  for  these  systems  "point  source  in  the 
crust"  will  be  defined  there. 

Section  4;  •  Preliminary  numerical  examples  for  the  case  of 
the  explosive  source  in  the  crust  will  be  pre¬ 
sented. 

2.0  The  displacement  potential  due  to  a  vertical  and  a 

horizontal  point  force  In  a  layered  medium  over  half- 
space . 

Throughout  this  paper  we  prefer  to  work  In  the  fre¬ 
quency  domain.  The  description  of  wave  propagation  by  a 
frequency  transfer  function  is.  In  many  cases,  more  advan¬ 
tageous  than  the  equivalent  description  by  the  system  re¬ 
sponse  to  a  unit  pulse  excitation,  since  In  man;'^  studies 
the  knowledge  of  the  frequency  response  Is  required  In  any 
event,  and  since  the  overall  transfer  functions  of  a  series 
of  linear  systems  can  more  easily  be  obtained  by  multi¬ 
plication  of  the  corresT>onding  transfer  functions  than  by 
convolution  of  the  corresponding  unit  pulse  responses. 

2 . 1  The  displacement  potential  due  to  a.  vertical  and  a 
horizontal  point  forcei 


We  shall  consider  the  process  of  wave  propagation 


6 


from  a  vertical  and  a  horizontal  point  force  in  a  layered 
medium  over  a  half-space  after  transforming  our  model  of 
the  process  iti  a  timeless  space  by  a  Fourier  integral . trans¬ 
formation  of  the  equations  of  motion  for  the  various  homo¬ 
geneous  layers  of  oui'  model,  including  the  source  and 
boundary  conditions.  In  such  a  space  information  from 
source  to  receiver  is  no  longer  transmitted  in  time.  In 
this  timeless  space  a  local  disturbance  of  tht  equilibrium 
at  the  source  will  "immediately”  affect  the  whole  system. 

Our  transformed  point  source,  a  Dirac -delta -source  in  space, 
which  acts  either  as  a  vertical  or  horizontal  point  force 
will  generate  a  disturbance  in  the  distribution  of  the 
stress-tensor.  The  stress-tensor  is  coupled  to  the  strain 
tensor  by  the  elasticity  tensor.  Therefore  the  redistribu¬ 
tion  of  stresses  will  also  cause  a  redistribution  of  the 
transformed  displacement  vector  throughout  the  transformed 
space,  right  to  the  place  of  the  receiver. 

For  the  description  of  our  model  we  choose  the  two 
coordinate  systems  depicted  in  figure  2.  For  the  most  part 
the  cylindrical  system  (r  ,0,Z  )  will  be  used.  In  Sect. 

2.2  the  right-handed  xyz-systera  is  sometimes  more  advantage¬ 
ous  . 

It  will  now  be  useful  to  define  transformed  displace¬ 
ment  potentials  fu  tl  and  from  which  the  transformed 
displacements  vector  (  derived  (see 

Harkrider,  1964)  as  follows; 
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1  Zxi 
^  r  a® 

(i.|) 

Vi(r,©,r)  »  f 

.  1  .  ®Xi. 

r  "57" 

(1.1) 

a*  ft  .  s, 

1-  Ti 

(1-3) 

la  Ij***  1  ^ 

We  will  then  carry  these  potentials  from  the  source  to  a 
"receiver"' within  the  homogeneous  half-space.  Comparing  these 
with  the  potentials  which  would  have  arisen  if  the  same  type 
Qf  source  had  been  placed  in  an  infinite  medium  of  subcrustal 
material,  we  arrive  at  the  quotient  of  the  two  potentials 
which  may  ;'>e  regarded  as  the  transfer  functions  of  the  source 
crust  system. 

The  potentials  will  not  be  carried  directly  from  the 
source  to  receiver.  The  Thoms on-Haske 11  matrix  formulation 
uses  the  so-called  motion-stress  vector,  whose  components 
are  the  components  of  the  displacement  or  particle  velocity 
vector  and  of  tbe  stress  tensor.  This  motion-stress  vector 
is  carried  across  the  Interface  continuously,  and  through 
the  layer  from  top  to  bottom  by  a  matrix  transformation. 

Thus  the  continuation  of  the  motion-stress  vector  from  source 
to  receiver  can  be  expressed  as  a  series  of  nwitrlx  transform¬ 
ations.  Dunkln  (1965)  reformulated  the  Thoms on- Haskell 
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matrix  method  to  reduce  numerical  difficulties  at  higher 
frequencies.  Dunkin' s  refined  matrix  formulation  allows, 
so  to  speak,  an  increased  flow  of  information  from  source, 
to  receiver  through  the  "noisy  channel"  of  a  digital  com¬ 
puter  . 


This  paper  will  follow  the  matrix  formulation  as 
explicitly  developed  by  Hcrkrider  (1964).  To  avoid  repeti¬ 
tion  we  will  simply  sum  up  that  he  continues  two  motion- 
stress  vectors  and  through  the  model  from  source  to 


receiver.  In  the  1 


th 


layer  these  vectors  are: 


Ki  *  {  ^  (2.k] 


"  { 


(a.sj 


These  vectors,  which  are  dependent  on  only,  are  essen¬ 


tially  derived  by  a  series  of  integral  transformations  from 
the  Fourier  transformed  displacement  potentials.  They  are 
connected  with  the  transformed  potentials  9iW,V8t*U;(*) 


’■  -  . . pitlit'liil-o 
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*  r  ,  d*s-iU)  .  1*  4  Vii^i  1 

««iW=  ^ 1 -3?— * TI*"  ^  TT^J 


-  Ai  Vi**) 


U ») 


and 


iliiil  •Lit  Xi**) 

r  .  lu-  ^x*‘*) 
tud)  s  *r‘  ji2 

In'*  eq.  (2.6)  -  (2.11)  the  following  notations  are  used; 


<110 


c 

Ki 

Si 


horizontal  wave  number,  with 
c  =  horizontal  phase  velocity 

dllatatlonall  wave  number  with  ^  , 

y  flti  =  p-veloclty  ln>  ^^h 

shear  J  =  S- velocity  InJ 


t/h 

Lame’s  constants  in  the  1  layer 
Prom  the  Fourier  transformed  displace¬ 

ment  potentials  fi  f  M'l  |  %i  are  obtained  by: 

fi  <  r|  ©!*/»/  fli^)  ^  I 


tilr,  ®,2)  »y  >rn*)  ‘••#6 

0 

»4..i©44 


(i.ia) 


(MJ) 


(J.M) 
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Jj^are  the  Bessel  functions  of  order  where  /  =  Of  I  ±n 
case  of  vertical  and  horizontal  point  force  respectively. 

This  paper  will  only  be  concerned  with  the  P-waves 
radiating  into  the  half-space.  Therefore  we  will  now  only 
follow  the  motion-stress  vector  from  source  to  receiver, 

This  vector  is  continued  from  top  of  layer  I  to  the  bottom 
of  the  layer  k  by  the  matrix  relation: 

^  *  'r-  I 

The  elements  of  the  matrix  C^n.  ,  the  layer  matrix,  are: 

s  «  Ifi  tmVi 

(••,)„  =K)«  [ifi  1, ^1 
(••il.  •KIm  ■  '‘‘1 


I* 


(4 
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where; 


Ifi  =  ^  (  c‘)* 

■  A  »«j  <ii 

Q;  S  A  r^.  di 

di  '  2l  -  *i-. 

I  -i  nrii 

\  nr^ 

i-i  /Fir 


/»* 


for 

A  ^  A^j 

for 

A  *  A^^ 

for 

A>  A^. 

for 

A  <  A|n 

(2I7J 


(i.iSj 
(2.  W) 
(2.20| 


(2.2I| 


12.22) 


Since  the  mot ion- stress  vector  is  continuous  across  the  in¬ 
terfaces  it  can  be  continued  directly  into  the  next  layers 
until  the  receiver  point  is  reached.  Harkrlder  then  shows 
that  the  point  forces  can  be  represented  in  our  Z  -dependent 
model  space  by  introducing  a  discontinuity  of  the  motion- 
stress  vector  at  the  place  ZsD  of  the  source.  These  dis¬ 
continuous  motion-stress  vectors  are  in  case  of  a  vertical 
point  source: 


where 


{  I 


I A 

Zw 


(2.2S) 


(2.24) 
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and  for  a  horizontal  point  force: 

{  0,  0,  0,  }  (a  w) 

where  —  • 

«rr», .  I W 

LM  is  the  Fourier  transformed  point  force.  Now  the  source 
layer  l»S  is  divided  into  two  layers  5,  and  Sj  at  the 
level  of  the  source  depth  D  (see  figure  2  ).  The  layer 

will  be  referred  co  as  the  top  source  layer.  The  layer 

matrix  for  the  layer  5,  are  the  same  as  for  the  undivided 

Sf  * 

source  layer  S  only  nas  to  be  replaced  by  the 

depth  of  the  source  in  the  source  layer. 

Taking  care  of  the  prescribed  components  of  the  motion- 
stress  vector  at  the  free  surfaces  and  the  radiation  condition 
that  no  sources  are  at  infinity,  one  is  in  a  position  to  con¬ 
tinue  the  motion-stress  vector  Introduced  at  the  source  to 
the  receiver  position  in  the  infinite  half-space.  In  the 
half-space  n  the  transformed  potential  in  the  form; 


Aj 

(2.27) 

.-“'••'-'"I 

MJ 

(2.28) 

j  «  V  or  H 


satisfies  the  wave  equatioiis  and  the  radiation  conditions, 
provided  and  are  chosen  as  in  ^?q.  (2.21),  (2.22). 

In  that  case  Aj  and  Wj  are  connected  with  the  motion-stress 
vector  /T|^  in  the  half-space  by  the  matrix  relation: 
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=  E|i  • 

•'K  *K 


j.  Vff  H 


(i-Mj 


where: 


CT** 

* 


-Mg) 

0 


0  0 

0  0 

I  0 


(iw) 


The  iOllowlnc  matrix  relations  connect  the  potential  vector 
defined  in  eq.  (2.26)  with  the  motion-stress  vector  (0) 
at  the  free  surface  and  the  scarce  motion-stress  vector: 
a)  for  a  vertical  point  force: 

i,  0 

,  J.  +  r:' .  0 


0 


{a.ji) 


\  V 

b;  -for  a  horizontal  point  force  in  dlrcctlon©«0  (see  figure 

2  ):  (a.  )  risdiil  r  0  1 


c 

c 

0 

0 


4  R' 


(1.32j 


where 


J  S  e;  . a 


>  It  ’ 


|2J3) 


is  the 
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will  be  called  the  crustal  inatrix 
>^erse  of 


and  where  R 


«•»! 


•  ••• 

n 

will  be  termed  the 

top  source 

matrix,  while 

is 

the  top 

layer  matrix. 

Harkrlder  ( 

1964)  derives 

a  very  con 

venlent 

form  for  the  elements  of 

(Rfc,,  )iiH 

0“*  - 

(  1 31 

I**n)i3 

("iwIh 

“  (R^silsi 

)ix 

*  I**®*  )<n 

)ll 

(UH) 


(2.»S) 


In  case  of  the  vertical  point  source  we  are  only  concerned 
with  the  -Tiatrlx  eq.  (2.31)^  for  the  horizontal  point  force 
with  eq.  (2.32).  Both  matrix  equations  represent  each  a 
system  of  4  linear  algebraic  equations  foi  the  4  unKnowns 

or 

tively . 

Up  to  this  point  we  have  followed  Harkrider's  develop¬ 
ment.  Since  he  was  in+’-erested  in  the  surface  displacements, 
he  solved  the  systems  of  equations  for  i 

arriving  finally  at  the  displacements  due  to  surface  waves 
of  the  Rayleigh  type. 

Since  in  this  paper  we  are  interested  in  the  radiation 
of  P  waves  into  the  half-space  we  shall  solve  the  systems 
for  Ay  and  An  ,  corresponding  to  P  waves  from  the  vertical 
and  horizortai  point  xirce.  Given  these  quantities  we  then 
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derive  integral  representations  for  the  transformed  poten¬ 
tials  due  to  a  vertioa]  point  source: 


(!?l 


and  nn  due  to  a  horizontal  point  force: 

9.»  *  r«.s(-'J2lM  (^1"  (14 

Solving  (2ol)  forAy  results  in: 


(iU) 


I 


(lit) 


where; 

U  [l  Jiri*n)(  hx^Jn^ii/  *  (  j  (*^*1  p3l*Ai)(Jii  Wn^ujj 

V  [(Jjj  )( 1  *  Ji4  ^  )  1*^*^  ^  ^  i) 

are  the  elements  of  the  crustal  matrix  J  (eq.  (2.33)), 
and  U  and  V  (using  eq.  (2.33))  given  by: 


) 

n-^oi 


V  -  <f^+  «»(»«* 

In  case  of  a  vertical  point  source  (see  eq.  (2.23), (2.24) ) 


we  have 

Uy 


(HJ 


^t. 


i/ 


UM) 

(2M) 


MiJ 
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and  for  the  horizontal  point  force  eq.  (?.4l)  and  (2.42) 
degenerate  to  (see  eq.  (2.26),  (2.26)): 

^  ‘W*  ^ 


U-  <'S| 


n 


Considering  that  is  a  common  factoi-  in  the  quantities 

U  and  V  and  therefore  a  factor  of  and  Ah  we  shall  define 
for  later  use  the  quantities  Ay  and  Ah  by: 


^  5v 

*  4u 


H  ® 


K 


The  eq.  (2.36)  and  (2.37)  represent  the  complete  solutions 
for  the  dilatational  displacement  potential  v/jthin  the  half- 
space  due  to  a  vertical  or  a  horizontal  point  force  jn  a 
layered  crust.  In  the  ne::t  section  ^r.2  we  will  use  contour 
integration  in  the  complex  il -plane  to  separate  modal  solu¬ 
tions  of  the  surface  wave  type  from  the  contributions  of 
body  waves. 

2 . 2  Separation  of  modal  surface  wave  solutions  from  body 
waves 

The  procedure  for  the  separation  of  modal  surface 
v;ave  contributions  from  body  waves  will  be  described  only 
for  the  transformed  dilatational  potential  due  to  a 
vertical  point  source.  The  separation  in  case  of  the  hori¬ 
zontal  point  force  follows  the  same  line  of  reasoning. 

The  dilatational  potential  in  the  half-space 


msr-fi 
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due  to  a  vertical  point  source  Is  given  by  eq.  {2. 36), which 
vie  rewrite  using  eq.  (2.45): 


I*'.!* 


r<«)* 

\  Utl 


where:  a  Z  *• 

and  Is  the  same  as  defined  in  eq.  (2.2l).  To  evaluate 

the  Integral  in  (2.43)  we  use  contour  Integration  in  the  com¬ 
plex  -plane.  vVe  will  Introduce  the  following  abbreviation: 


U  i  /  A 


=  ^  CwT  Avui 


(IMI 
(t  WI  • 


IX.SC) 


The  Be.gsel  function  will  be  replaced  by  the  Hankel  functions: 

JolH  H.'‘'(Ar)] 


With  this  equation  (2.48)  ^nay  be  written  as; 


•  I,  ♦  ii 


where: 


I. 

It 


e  • 


di 

Ho^^^{kr)di 


Prom  now  on  we  consider  the  wave  number  m,  c'^d  the  angular 
frequency  CJ  to  be  complex  variables: 


t;  -- 

(0  «  s  -  1 1 


(ISIJ 

(ISi) 


IJSS) 

(»S*J 


I 


IS 


Keeping  in  mind  that 


lU  for  4j'0,t>0 

^409 

tu.  H.‘^|?r)^0  for  A>0,t<0 


(J.srj 

(i.SIJ 


v^e  displace  the  path  of  Integration  from  the  positive  real 
axis  to  the  imaginary  axis  for  the  integral  I|,  containing 

,  and  to  the  negative  imaginary  axis  for  the  Integral 

U 

containing  ,  closing  the  path  by  the  qua’^ter  circle  in 

the  first  and  fourth  quadrant,  respectively.  The  contribu¬ 
tion  to  the  integrals  X|  and  I2  along  these  quarter  circles 
vanish  according  to  eq.  (2.57)  and  (2.58). 

The  chosen  deformation  of  the  original  contour  is  Indi¬ 
cated  in  figure  3.  The  integrands  in  Ij  and  I|^  have  poles 
and  are  multivalued  functions  of  ^  ‘  The  poles  occur  at  the 

zeros  of  the  denominator  of  There  is  a  singularity 

at  ^*0  for  the  Hankel  functions.  This  singularity  dis¬ 
appears  as  both  Integrals  are  summed  up,  since  the  Bessel 
function  is  regular  at  It  will  therefore  be  disre¬ 

garded  . 

Due  to  the  presence  of  the  expressions  i 
v/hlch  appear  in  the  crustal  matrix  J  ,  the  integrand  is  mul¬ 
tivalued.  To  get  a  single  valued  function  we  must  intro¬ 
duce  2t<.  branchllncs  originating  at  the  branch  points 

(  ^  *  ""I  •  To  assume  the  vanishing 

of  the  potential  for  we  must  remain  on  that  sheet  of 
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thf^  Riemann  surface  where 

1U(i  )>0  ;  >0 

t  s  l|  • .  • .  |K 

Therefore  the  branch  lines  are  given  by: 

V(i^rpjJzO 

t  I  K 

The  positions  of  the  poles  and  branchlines  are  Indicated  in 
figure  3  for  the  case  of  lli(<o)>0  .  For  the  poles 

and  branchlines  would  be  displaced  to  the  first  quadrant.  But 
the  results  of  the  following  derivations  are  the  same  in  both 
cases. 

Jardetzky  (1953)  pointed  out  that  in  all  cases  of  wave 
propagation  from  a  point  source  in  a  multilayered  medium  over 
the  half-space,  all  expected  branch  line  integrals  vanish 
except  one  branchline  in  each  potential  ^  and  y  correspond¬ 
ing  to  the  terras  or  , respectively,  for  the  half-space. 
(See  also  Ewing  et  al.  (1957)).  This  can  also  be  proven  by 


(l.SSj 


considering  the  symmetry  properties  of  the  layer  matrices 

(see  eq.  (2.16))  and  Ea  *  (see  eq.  (2.30)).  In  this  way 
it  is  possible  to  prove  the  vanishing  of  the  branchline  inte¬ 
grals  due  to  I  for  i  s  |  .  There  still 

remain  two  branchline  contributions 

and  of  the  half-space  for  both  dJ latational  and  shear 
potential . 


(Z(>0) 

(Ml) 
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Referring  to  figure  3  the  contributions  to  the  Integrals 
I|  and  from  branchlines  and  poles  may  be  written  as  follows 

»  i  fjiiz)  H,‘"  i;Tf» 

o  ^ 

<•*  L  IJ  I  ▲ 

n®  (  I  W  rj  «  t  T 


i/...,  -i/....  -  iiri  Z 

<‘K 


(iti) 


lb  can  be  shown  that 


Ti  iT)«-  T(-ir| 


Considering  further  that 

II  t J  *  « _ 


12.4SJ 


we  may  write  X  |  as: 

I,  «.'*'(  itr^T  ( 

(I 

Summing  1|  and  Ij  according  to  eq.  (2.52)  we  realize  that  the 
contributions  along  the  positive  and  the  negative  imaginary 
axes  cancel.  Therefore  the  final  complete  solution  for 
is: 

Up  to  eq.  (2.6?)  no  approximation  has  been  used.  Thus  eq.{2.67) 


(2W| 


(1.47) 
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represents  the  exact  solution  for  the  r'ourler  transformed 
dilatattonal  potential  in  the  half-space  due  to 

a  vertical  point  source  within  the  crust. 

We  can  now  separate  body  v/aves  from  surf  a  e  wave  con¬ 
tributions  by  the  following  line  of  reasoning.  Body  waves 
may  be  regaraed  as  those  waves  which  travel  in  the  half¬ 
space  with  horizontal  phese  velocity  .  All  waves 

traveling  with  phase  velocity  C5|\ii  are  more  or  less 
trapped  in  the  crust,  l.e.  the  amplitude  diminishes  rapidly 
within  the  half-space  with  increasing  depth  Z  .  They  are 
mathematically  described  by  the  residue  contributions  from 
the  poles.  Thus  v;e  may  regard  the  contributions  from  the 
branchline  contours  as  descrloing  the  body  waves  in  the  half- 
space  generated  by  the  point  source  in  the  layered  medium. 

We  now  define  the  transformed  dilatatlonal  potential  of 
body  waves  vy  in  the  halfspace  due  to  a  vertical  point 
source  by:  ,  ^  - 

This  definition  of  body  waves  is  especially  useful  at  large 
distances  from  the  source  such  that  the  suri  -  oe  wa'^^es  have 
had  sufficient  time  to  separate  from  the  boay  waves.  We 
shall  now  restrict  our  attention  to  the  value  of  the  Inte¬ 
gral  representation  of  at  large  distances  from 


(i.ig) 


i 


the  source. 


pp 


2 . 3  Large  distance  approximation. 

At  iarr:e  distances  iron  the  source  v/c-  substitute  the 
approximation  of  the  Hankel  Vinctior.s  for  larje  vaiues  of 
the  argument  into  eq.  (2,1),  ’..■riting  the  two  integrals 
in  one  as: 


- ,  ,  f  "'I 


Por  large  distances  v/e  shall  approximate  the  value  of  this 
Integral  by  the  contribution  at  the  saddle  point.  To  apply 
the  saddle  pojnt  method  we  revirlte  the  integral  in  the  form: 


Vv’here 


f(S  ) «-‘  ( ^  1) 


(i.70) 


(1.71) 


Since  we  are  considering  only  the  branchline  Integrals  where 
there  are  no  poles  of  m)  on  the  chosen  branch  of  the 
Riemann  surface,  we  are  Justified  in  assuming  that  the  factor 
of  <«f  m]  varies  slowly  compared  to  itself, 

especially  for  large  f  .  In  the  fourth  quadrant  the  saddle 
point  of  the  Ir.tegrand  is  determined  from  f  to  be 

located  at; 


•wnere 


«  =  'I  ♦ 


(jw) 


I2.7ij 


Is  the  distance  from  the  plurr/o  point  uf  the  source  at  tne 


base  of  the  crust  to  the  receiver  at  (VjZ)  .  Since 


always 


one  may  conclude  that,  as  long  as  K  is  sulTlciently  greater 


than  f  ,  l.e.  as  long  as  the  receiver  is  suificiently  re¬ 


moved  from  the  base  of  the  crust,  the  saddle  point  is 


separated  from  the  branch  point  at  • 


Writing  he  complex  branchpoint  Stt|^  as: 

*  K-k  I 


the  path  of  steepest  descent  tlirough  the  saddle  point,  on 


which 


has  a  maximum  and  decreases  most  rapidly 


has  the  direction  ^  ; 


♦  =  f  •? 


Evaluation  of  (2.69)  at  the  saddle  point  along  the  path  of 


steepest  descent  yields  finally 


)  r  T(l,)  i  e 


Here  we  ha^o  returned  to  real  angular  frequency  W  and  real 
wave  number  A  .  Slc.ce 


X  «  m  ^ 


and  since 


I  a 

*«(h  =  «  ■ 


the  saddle  point  is  located  at: 

4,  =  ^  ^  1  = 


(iUH) 


IJ.75) 


(2.77) 


(2.7»j 


(2.7^) 


f2t0) 


Together  v;itn  eq.  (2.30)  -'or  T(ll)  vje  .nay  thiarei’ore  vnr'lte: 

_  .  .14  ..  A 

{—I  —J - 


v;here  jse  has  been  made  oi': 

I  - 

is  t’ne  angle  of  incidence  at  v;nich  the  P-v;ave3  leave  the 
base  of  the  crust  (see  Tigure  l). 

The  equivalent  result  ior  the  dilatational  potential  due 
to  a  horizontal  point  i'orce  v;ith  direction  is: 

-1^1  •  AgCayO  — 

<p»,»(r,©,ij=  ty  Iw)  n 


In  the  next  section  we  shall  naw  derive  the  expressions  for 
three  types  of  point  forces  commonly  used  as  models  in  focal 
mechanism  studies,  namely:  an  explosive  point  source,  a  single 
couple  and  a  double  couple. 

3.0  The  transfer  function  for  P-waves  at  large  distances 
Trom  thjee  types  o7  point  sources:  explosive  source, 
single  couple,  and  double  couple. 

3 . 1  Procedure  and  general  remarks. 

In  this  section  we  shall  derive  the  transformed  po¬ 
tentials  at  large  distances  for  the  following  types  of  point 
sources: 


1. 

Explosive  point  source,  j 

potential 

2. 

Single 

couple,  potential 

in  3.3 

3. 

Double 

couple,  potential 

in  3.^ 

(Ml) 

(MJj 

|2.«) 


These  multipole  types  of  point  sources  can  be  derived  from 
a  single  point  force  by  spatial  differentiation  and  super¬ 
position  (Ewing,  et  al.  1937).  We  shall  take  advantage  of 
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the  fact  that  the  displacement  field  due  to  a  point  force 
of  arbitrary  orientation  in  a  layered  system  can  be  com¬ 
posed  by  superposition  of  the  fields  due  to  the  normal 
(vertical)  and  the  tangential  (horizontal)  component  of 
that  point  force  (Keylis-Borok,  1953)- 

We  wish  to  c'^mpare  the  potentials  of  the  point  source 
in  the  layered  medium  and  of  the  point  source  in  the  in¬ 
finite  medium  of  subcrustai  material.  The  latter  potential 
we  shall  denote  by  ,  where  i  -  *I^|3  depending  on 

the  source-  type.  We  shall  define  functions: 

,  4L 

which  we  shall  call  transfer  functions.  They  can  be  inter¬ 
preted  as  describing  the  effect  of  the  following  operations 
on  the  particular  type  of  point  source  in  the  infinite 
homogeneous  medium  of  subcrustai  material: 

1)  Replacing  the  half-space  on  top  of  the  point 
source  by  a  system  of  (n-l)  different  homogeneous 
layers  with  paralell  interfaces. 

2)  Moving  the  point  source  into  the  s^'^  layer,  the 
source  layer. 

Therefore  we  shall  call  for  t  : 

T^  :  dilatational  transfer  function  of  the  system 

f 

"type  t  point  source  in  the  crust"  at  large 
distances . 

TjL  has  been  defined  as  the  transfer  function  for  the  dilata¬ 
tional  potential.  It  should  be  mentioned  that  any  quantity 
derived  from  the  dilatational  potential,  such  as  displacement. 


. . Ilium . . . . . . .  ■"  '  . . .  . . .  . . . 
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particle  velocity,  etc.,  has  the  same  transfer  function 
as  long  as  one  is  operating  in  the  same  way  on  the  'input" 
and  the  "output 

For  all  three  types  of  point  sources  cr">sidered  here 
the  following  explanations  apply: 

1)  A  point  force  may  act  in  the  arbitrary  direc¬ 
tion  of  the  unit  vector  K  ,  which  may  be  defined 
by  its  components  in  our  xyz-coordinate  system: 

^  •  {  *^1  4  I  •'a  J 

We-  shall  decompose  the  point  force  ^  into  its  com¬ 
ponent  vectors: 

Ki  *  •(  »|  K I  ®  I  *  ) 

{  0 ,  Hi  K  I  0  ) 

i?,  =(0,  0,  n,  Kj 


I VI) 


(V)| 

a- 4) 

(IS) 


where  K  ®  I  ^  I 

We  are  now  dealing  with  one  vertical,  ,  and  two 
horizontal,  point  forces,  and  .  The  potentials 
and  fiH  due  to  these  vertical  and  horizontal 
corrponents  are  obtained  from  the  expressions  (2.8l) 
and  (2.83)  derived  in  section  2.3  by  taking  into 
account  the  decompositicii  of  the  force.  This  means 
that  in  the  motion-stress  vector  (eq.  (2.23)  and 
(2.25))  we  have  to  replace  ^  by  ^  or  or 
^**1  for  the  thrae  components  of  the  point  force. 

2)  The  potentials  and  fH  due  to  decomposed  couples 
of  forces  (that  is,  a  pair  of  equal  but  opposite 
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forces  acting  at  che  ends  of  an  Inf Initeslmal 
line  element  ,  where  is  the  unit  vector  in  the 
direction  of  the  line),  these  potentials  are  obtained 
from  and  oy  taking  the  component  of  the  grad¬ 


ient  of  these  potentials  in  the  direction  of  (see 


lor  Instance  Nakano  (1923)): 

y'l)  a  tiiS  •  •J  V»  f  H 


(V7i 


The  subscript  affixed  to  the  gradient  symbol  indi¬ 
cates  differentiation  with  respect  to  the  source 
coordinates.  The  distance  from  the  plumb  point 
of  the  source  at  the  top  of  the  mantle  to  the  re¬ 
ceiver  point  as  defined  in  eq,  (2,73)  does  not  change 
on  vertical  variation  of  the  source  location.  There¬ 
fore  the  derivative  f  R  with  respect  to  the  vertical 
source  coordinate  vanishes. 

3)  The  gradients  in  (3.6)  and  (3.7^  require  operation 
on  Av<5)  and  Adh)  defined  in  eq.  (2.46)  and 
(2.47).  It  should  be  realized  that  these  two  quanti¬ 
ties  are  only  dependent  on  the  vertical  of  the  source 
coordinates,  while  derivatives  with  respect  to  the 
horizontal  source  coordinates  will  vanish.  The 
crustal  matrix  J  (see  eq.  (2.33))  does  not  contain 
the  source  depth  2>  as  parameter;  only  the 

expressions  forU  andV  are  t  dependent,  since  they 
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contain  the  top  source  maLrix  (see  eq.  2.3-4) ); 

Br..  =  . 


'Rj,  =  -H  "S-l  .  "I 

In  this  matrix  product  only  the  top  source  layer 

matrix  0|,  is  ^  dependent.  We  define  tne  deriva¬ 
tive  matrix  ft'jj  ; 

and  likewise; 


*  ••••  • 


or 


®  Hsi  • 


Referring  to  eq  (2.4l)  and  (2.42)  for  U  and  V  we  need 
only  the  following  elements  of  : 


The  elements  of  0|,  which  will  be  used  for  the  com¬ 
putation  of  these  4  elements  are  (compare  eq.  (2.l6)): 


(•tllll  * 

fh 

r,,  aim  f, 

>  -(h- 

1 )  rp.  *i« 

1.1  = 

U  [(ft 

•  1  j  Ci>  Pj  ' 

♦  1*  •‘p) 

e««,] 

K.1 

= 

r^,  ail. 

(•«: 

^ 

on  Qj) 

(2.,1 

111  " 

1  «>  tti3 

111  - 

^  [di*' 

)r^  1^1 

(•'»')  ij  * 

u>  0*] 

)aii » 

(311 


f  v<  I 


(3I0I 

Cl-irl 


(in) 


/r0'' 
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The  quantities  ^  ^ 

Sv  *|^4v_  'al 

are  obtained  from  Ay  and  A||  by  simply  replacing  the 
elements  of  the  top  source  layer  matrix  ft||  by  those 
of  the  derivative  matrix  A*»i  . 

4)  For  handy  reference  in  the  following  sections  we  de¬ 
fine  the  following  two  vector  quantities: 

^  ^  cmO  &v  ;  ^  I  ®  *'  } 


|V 


(J.litl 


If  .  Uc*^*;**®  a'hJ  U  'S) 


5)  The  transformed  dilatational  potential  due  to  a 
point  force  acting  in  direction  ft  in  an  infinite 
homogeneous  isotropic  medium  is  given  by: 


(Ml) 


where  is  the  density  of  the  infinite  medium  made 
up  of  subcrustal  material.  At  large  distances  this 
is:  ^  ^ 

f--  T— 

Here  R®  is  the  unit  vector  in  the  direction  from 
source  to  receiver: 


(JI7) 


llldWliii 
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From  we  derive  the  potential  due  to  a 

couple  of  forces  acting  on  a  line  element  with  direc¬ 
tion  by  determining  the  comnonent  of  the  gradient 
of  in  the  direction  (compare  eq.  (3*6)  and 

(3.7)): 


*  •* 


(Jl'i 


This  is  for  large  distances:  ^  --‘Ik* 


fU  » 


3.2  The  explosive  point  source. 

An  explosive  point  source  can  be  represented  by  3 
mutual  perpendicular  couples  of  equal  point  forces  without 
moment.  In  case  of  the  layered  medium  we  will  represent  the 
source  by  one  vertical  couple  in  the  vertical  z-direction 
and  two  hoi*izontal  couples  in  the  x-  and  y-directlon.  The 
potentials  for  these  couples  are  obtained  by  operating  with 
(3.6)  and  (3.7)  on  the  expressions  for  and  fn  from 
eq.  (2.8l)  and  (2.83).  For  the  dllatational  potential  of 
the  vertical  couple  we  obtain  at  large  distances,  omitting 
terms  of  higher  order:  ... 

»'»  (r,7|  R 

likewise  for  the  sum  of  the  potentials  due  tc  the  two  hori¬ 
zontal  couples:  ^ 

v'h  «  ‘  (§)  “r 

•  I  <-.) 

Superposing  the  two  potentials  fyand  f ^  we  finally  arrive 

at  the  dllatational  potential  due  to  an  explosive  point 


(S.22) 
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source  in  a  layered  medium: 


The  same  source  in  the  infinite  mediun:  of  subcrustal  mater¬ 
ial  gives  rise  to  the  following  potential  at  large  dis¬ 
tances: 


Vi- 


-I 


This  may  be  realized  from  (3.20)  by  remembering  that  for  all 
three  couples  without  moment  l\  is  equal  to  ^  .  Thus  essen¬ 
tially  the  following  terms  have  to  be  summed: 

? .  ?  I  ^ 


( +(ir,-  R*)* 
+  in)S  («•,)"  = ' 


(Vix) 


(JiSj 


According  to  (3.1)  the  transfer  function  T^^  is  then 

T,  »  2*1^^  j*  [- j  i'y  Am  ] 

where  use  has  been  made  of 


3.3  The  single  couple  type  of  point  source 

The  single  couple  type  of  point  source  is  depicted  in 
figure  4.  Without  loss  of  generality  we  shall  choose  the 
X  I  Yi  Z *  system  orientation  so  that  the  horizontal  component 
of  the  force  is  in  the  direction  of  the  x-axls  (  ©*  • 

The  direction  of  unit  vector  Iah  along  the  line  dement  is 


arbitrary,  the  only  constraint  bein'  that  the  point  force 
end  the  line  element  are  perpendlculart 

55  •  if  «o 

The  following  notations  will  be  used: 

(p  ;  plunge  of 

^  •  trend  of  ^ 

i  plunge  of 


(Ml) 


Thus  and  have  the  following  components  in  the  xyz-system: 


,  0,  J 

The  cons'traint  (3 *26)  gives  rise  to  the  following  relation: 

4.  viklf  sO 

We  shall  distinguish  two  cases; 

Oase  A:  <jp^0®or  IfO^  and  or  270® 

Case  B;  <|paO*  or  ltd*  and  or  110^ 

For  case  A  relation  (^.31)  is  satisfied  if: 


(J.iO) 


{j»i) 


-- 


~  - 


««»<  f»X 


(lU) 


(1  Jl) 


. 


\vnere 


iJ-  *  .  V  •  -  vU'x  t*»*V  +0 


In  this  case  the  two  unit  vectors  M  and  have  the  lollow- 


ing  components: 


» 


{vi*yt»>X,  ‘^X  I 


T?  r  <  ,  0,  ti»f  J 


(l.Jf) 


(134) 


For  case  B  relation  (3-31)  iJ^  satisfied  for  arbitrary  plunge 
of  ivT  .  The  two  unit  vectors  then  become: 

iSf  «  <  0  ,  mj-A,  Ui-A  J 

*v  «  {  •  I  ®  I  0  ] 


t3.17) 

IJ3S) 


Referring  to  the  remarks  on  the  decomposition  of  force  couples 
in  section  3.1  and  applying  eq.  (3.h)  and  (3.7)  to  the  poten¬ 
tial  of  the  vertical  and  horizontal  component  forces,  we  find 

-I  I 

for  the  dilatatlonal  potentials  fztt  large  dis¬ 

tances. 

Here  V  and  H  are  the  vectors  defined  in  (3*1^)  and  (3.15). 
The  sum  of  both  potentials  is  the  dilatatlonal  potential 
of  a  point  source  of  the  single  couple,  t^^-pe  in  a  layered 


(3)1) 


(310) 


rnedlum; 


The  same  source  in  the  infinite  medium  of  subcrustal  mater¬ 


ial  produces  the  following  dilatational  potential 
(see  eq.  (3.20)): 

From  eq.  (3*1)  the  transfer  function  To  is: 

T  .  ?  4rii,  [-iiMiiyiig-y  )♦<«)»  «>©(«»•?)] 

3.4  The  double  couple  type  of  point  source 

The  so-called  double  couple  type  of  point  source  is 
equivalent  to  a  pair  of  tension  and  pressure  couples  without 
moment  (see  e.g.  Stauder,  1962).  We  shall  denote  the  unit 
vector  in  the  direction  of  tension  by  i  ,  and  the  unit  vector 
In  the  direction  of  pressure  by  p  .  Both  are  mutually  per¬ 


pendicular; 


t  -p  s 0 


Without  loss  of  generality  we  will  choose  our  x»yi2  system 
so  that  the  tension  couple  is  always  in  the  XjZ  -plane.  The 
plunge  of  t  is  denoted  by  y  ,  therefore: 

t  ■<‘»f ,  ®  I  } 

The  orientation  of  H  will  be  described  by  its  plunge  V  and 
its  trend  ^  : 


UUl) 


(MM 


(J.45| 


b<*4) 


The  con 

St  '1  ’.nt 

(3. 

T4)  gives  r 

is*- 

to  th 

e  sar 

ne  relation  as 

in  eq. 

(3,31). 

A  s 

in  the  pre 

vie 

.us  cas 

a  vju 

must  distinquish 

case  A 

and  case 

p, 

no* 

Case  A: 

«p  +0* 

or  ♦  110*, 

X 

+40* 

or 

?  * 

4_  ;vi.f  t*»i( 

iUt 

,-«»♦  c*>x  ^ 

(V47I 

?  . 

( 

U>.f ,  0 , 

!«.»  ) 

IV4JI 

where  sUi  4  is  the  same  as  in  eq.  (3-3^) 
Case  B;  sO®  or  f  or  ^*^0® 

{  0, 

t  .(•,  0.0} 


(^.41) 

(i-so) 


By  the  same  method  as  In  the  previous  section  we  derive  the 
potentials  due  to  the  vertical  and  the  horizontal  componento 
of  the  tension  couple  from:  i  lA  |1 

viv  *  (^?r " 

;  _  I  till 


For  large  "R  »  I  v/e  neglect  highe-  order  terms  of  the 
gradient  .. 

Vtv  *  M  Tu  )  ^  I  ^ 

( lift fM  Ctiif  — “  un  0  ( ^  ^ )  s 
ttii  «  •  \  w /  K 

V  and  M  are  again  defined  in  eq.  (3.1^)  and  (3-15). 


(V5I| 

(V52) 


In  the  same  manner  we  derive  the  potential  due  to  the  vertical 
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component  of  the  pressu.’e  couple: 

(s)‘  «*»  ‘ . 

where  p»  Is  the  7-Gomponent  of  ^  ,  eq.  (3.47)  or  (3.49). 

The  two  potentials  due  to  the  x-  and  y-component  of  the  pres¬ 
sure  couple  are: 


fi.  ®  (irr  '•«.  P«  ^  n 


IVSS) 


(l.Sfc) 


(VSfJ 


Both  potentials  can  be  combined.  This  gives 


in  case  A; 


and  the  potential  is  in  this  case: 


«•*  •  f  •In  1^  ^ 

?|IV  *  *  Ul  )  ** 


tm 


9fimt  14 


IT 


if-^) 


|VS«) 


|3,S<) 


and  in  case  B: 

I 


l4  /-*  tTi 


e-iM 


(3'fcO) 

(iU| 


The  resultant  dilatatlonal  potent j  al  ^  of  the  double  couple 
is  the  sujn  of  four  potentials: 


4^  4 


.  I 


fi  -  V«v  ♦  tjH  »p»  ■*  »|ilt 


(l.Uf 
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We  distinguish  again  the  two  cases  A  and  B. 

Case  A:  With  eq.  (3.53)j  (3.34),  (3.5&),  and  (3.39)  we  de¬ 
rive: 


27  [ 


♦  S?2u. 


(x-eKf-tTil 


(Vfc») 


Case  B:  Prom  eq.  (3.h3),  (3.5^),  (3.60)  and  (3.61) 

? j e.z ) .  ( %  [-L  -u.  © ( r-  r ) 


The  same  double  couple  source  In  the  infinite  medium  of  sub- 
crustal  material  produces  the  following  dilatatlonal  poten¬ 
tial  • 

‘IB 

(3.65)  Is  derived  according  to  (3-20)  by  summing  the  effect 
of  the  tension  and  of  the  pressure  couple.  The  transfer 
function  T^  la  derived  with  (3.1).  Distinguishing  the  two 
cases  we  find: 

Case  A; 


T,  -  j> 


(VfcSI 


(J-u) 


!•(««** f-t  f'V  )-C»»fl(  t  M  )4io>^<^9(  f-  if  ij 

w*  (?•?}*  . 


3.3  Remark  on  nodal  planes 

The  single  and  double  couple  type  of  point  source  in 
the  infinite  medium  have  nodal  planes  where  the  potential 
vanish.  These  nodal  planes  are  defined  by: 

a)  Single  couple  (see  eq.  (3.42)): 

)(  sr  s^i  *0  (J  ‘«) 

that  is 

a  0  1S  •  R*  »0  (5’^^  I 

The  potentials  are  zero  in  planes  normal  to  “J? 

(the  direction  of  the  force)  and  normal  to  wt 
(the  direction  of  the  line  element). 

b)  Double  couple  (see  "q.  (3.65)): 


(?'.?•)*  s  ( 

u-yof 

This  is  fulfilled  If: 

p  i*  «  * 

( i-yi  1 

or 

0 

0 

Since  ^  and  t  are  unit  vectors  their  sum  and 
their  difference  form  vectors  in  the  direction 
of  bisecting  lines  betv/een  the  pressure  and  the 
tension  couple.  (3.72)  states  that  the  poten¬ 
tial  vanishes  in  the  planes  ixormal  to  these  two 


bisecting  lines. 


^9 


On  these  nodal  planes  the  present  definition  of  the  trans¬ 
fer  function  T2  in  eq.  {3. "^3)  ard  (3.66)  and  (3 .67) 

lose  their  meaning,  since  here  the  denominator  becomes  zero. 

In  general  it  cannot  fc  expected  that  the  dal  planes  of 
the  source  in  the  infinite  medium  do  coincide  with  the  nodal 
planes  for  the  point  source  in  the  crust.  .-Je  will  by-pass 
this  problem  by  defining  modified  transfer  functions  and  Tj  { 

ft  « Ti 

and 

These  modified  transfer  functions  will  not  "blow  up"  on  the 
nodal  planes. 

The  explosive  point  source  in  the  infinite  medium  does 
not  display  any  nodal  lines.  Therefore  there  is  no  neces¬ 
sity  of  defining  a  modified  transfer  function  fn  that  case. 


4.0  Preliminary  numerical  computations  for  the  transfer 
function  of  an  explosive  point  soured 

4.1  Introduction 

mm  0^ 

The  transfer  functions  T|^  and  1^  (or  Tj  and  Tj  ) 
are  given  in  a  form  which  is  convenient  for  programming  for 
automatic  c on^ju tat ions .  The  three  functions  have  been  pro¬ 
grammed  in  IBM  Fortran  II  for  the  IBM  7072.  The  transfer 
functions  contain  not  only  the  parameters  of  the  source  in 
the  infinite  medium  but  also; 

a)  the  parameters  of  the  crustal  model: 

I  Pi  I  Jl  f  ’ 


t  •  « 


I 


•  i^the  layer  thicknessesj 

b)  The  depth  ^  o!'  the  source  in  the  source  layer 

c)  The  angle  of  incidence  ^  into  the  mantle. 

In  this  paper  vje  will  only  give  preliminary  numerical 
results  for  the  explosive  point  source  as  the  simplest  t;/pe. 
There  is  no  azimuthal  dependence;  no  nodal  planes  exist  fOx- 
the  source  in  the  Infinite  medium.  The  selected  examples 
are  not  considered  to  be  complete^  but  they  will  serve  as 
guide  to  further  detailed  computation. 

liven  in  the  case  of  the  explosive  source  we  will  confine 
ourselves  to  the  study  of  the  following  questions: 

1)  Given  a  standard  crustal  models,  how  does  the  source 
depth  Influence  the  transfer  function? 

2)  For  a  fixed  source  depth  in  the  standard  model, 
how  does  the  transfer  function  vary  with  the  angle 
of  incidence  ? 

3)  Hov;  do  minor  variations  of  the  standard  model  In¬ 
fluence  the  transfer  function? 

The  program  for  the  transfer  function  Tj  is  based  on 
eq,  (3.26)  In  section  3*2.  For  the  computation  of  Av 
and  Ah  the  elements  of  the  crusta  .  matrix  J  .  (2.33), 
and  the  quantities  U  and  V,  eq.(2.4l)  and  (2.^8)  must  be 
determined.  With  some  minor  modifications  the  program  uses 
the  same  techniques  of  multiplication  of  matrices  with  real 
and  imaginary  elements  as  thiat  of  Hannon  (1964a, b). 
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The  program  has  been  tested  in  the  case  of  an  explosive 
point  source  by  comparing  the  theoretically  known  response  to 
a  4 -excitation  for  the  point  source  in  a  homogeneous  half- 
space  with  the  response,  calculated  from  the  source-crust 
transfer  function  T  ,  by  Fourier  synthesis,  of  a  layered 
model  in  which  all  layers  have  ldentlc'^1  properties.  For  the 
models  used  in  this  section  the  program  has  bee'^  r  n  for  a 
frequency  range  0  to  5  cps  without  "blowing  up." 

4 . 2  The  Influence  of  the  angle  of  incidence  y 

We  will  place  the  explosive  point  source  at  a  fixed 
depth  of  30  km  in  a  standard  crustal  model,  the  same  as  used 
by  Van  Nostrand  (1964): 

Table  1:  Standard  Crust  NOCR 


Depth  interval 
(km) 

P- 

( km/sec ) 

S-Velocity 

(km/sec) 

Density 

(g/cm3) 

0-3 

3.10 

1.79 

2.35 

3-20 

6.03 

3.48 

2.60 

20  -  45 

6.97 

4.02 

3.00 

45  -  00 

8.10 

4.67 

3.35 

The  crustal  model  is  displayed  on  the  right-hand  side 
of  figure  5.  The  amplitudes  of  the  transfer  functions  have 
been  determined  for  the  following  angles  of  Incidence: 

J  s  IS®  I  50*  j  45®  I  bO®  I  75®  ,  The  values  of  the  corres¬ 
ponding  horizontal  phase  velocities  C  are  indicated  in  the 
figure.  It  may  be  concluded  that  t'-a  spectra  of  P-waves 
incident  into  the  mantle  at  various  angles  do  not  vary  sig- 
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nlflcantly  with  angle  of  incidence  in  a  frequency  range  from  C 
to  about  0.3  cps .  This  is  especially  true  for  angles  of  inci¬ 
dence  between  15°  and  45^.  The  first  peak  at  O.O5  cps  and 
the  first  trough  at  0.1  cps  can  be  attributed  to  constructive 
and  destructive  interference  between  the  direct  P-wave  and 
the  pP  reflection  from  the  free  surface.  In  this  part  of  the 
spectrum  the  transfer  function  is  predominantly  governed  by 
the  Interference  between  P  and  pP  waves.  The  fine  struc¬ 
ture  of  the  crust  has  practically  no  influence  on  the  spec¬ 
trum.  The.  slight  shift  of  peaks  and  troughs  towards  higher 
frequencies  at  increasing  angle  of  incidence  may  be  accounted 
for  by  decreasing  time  lags  of  reflections  within  the  crust. 

4.3  The  Influence  of  the  source  depth 

Using  the  same  crustal  model  as  in  figure  5  we  will 
now  place  the  source  at  three  different  depths:  0.5,  8.0  and 
30.0  km.  The  angle  of  incidence  into  the  mantle  is  chosen  as 
^  =  45°.  The  amplitude  of  the  transfer  functions  are  dis¬ 
played  in  figure  6  for  the  frequency  range  0  to  O.3  cps.  For 
a  source  depth  of  0.5  km  the  first  peak  due  to  constructive 
P  and  pP  interference  should  occur  at  1.6  cps.  Therefore 
I  T|  is  relatively  flat  up  to  about  0.2  cps.  This  becomes 
clear  on  comparison  with  ITI  for  O.O  km  depth  and  30  km  depth 
where  the  first  peaks  are  located  at  0.2  cps  and  O.07  cps, 
respectively. 

It  may  be  concluded  that  the  transfer  function  for  an 
explosive  source  in  a  crust  is  rather  sensitive  to  changes 
'in  the  depth  of  the  source.  This  becomes  even  more  evident 
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J  if  one  considers  the  spectrum  up  to  4  cps. 

In  figures  7  and  3  the  amplitude  |T|  and  the  phase  ^  of 
is  depicted  for  two  source  depths:  Oo  and  S.O  km  for  J  =  45° 
Note  especially  the  very  high  amplicudes  for  the  near  surface 
source  between  1.2  and  2.0  cps,  and  the  relatively  low  ampli¬ 
tudes  fur  very  low  frequencies  between  0  and  0.2  cps  and  the 
band  between  3.0  and  3*4  cps.  Comparing  this  with  the  spec¬ 
trum  of  the  deeper  source  one  may  draw  the  following  tenta¬ 
tive  conclusions: 

As  the  source  moves  closer  to  the  free  surface  the  low 
frequency  part  of  the  transfer  function  becomes  a  band  of 
relatively  low  amplitudes  compared  to  the  amplitudes  of  the 
higher  frequency  bands;  e.g.  for  a  point  source  at  0.5  ion 
depth  in  the  standard  model  NOCR,  the  average  amplitude  be¬ 
tween  0.0  and  0.3  cps  is  about  1/5  of  the  average  amplitude 
in  the  interval  1.2  and  1.5  cps.  In  this  case  the  system 
acts  effectively  as  a  rejection  filter  for  the  lower  fre¬ 
quencies.  The  closer  to  the  surface  the  source  is  located  - 
without  violating  the  free  surface  condition  -  the  broader 
is  this  rejection  band.  This  range  of  small  amplitudes  is 
rather  flat  con5)ared  to  the  spectrum  for  the  same  range  of 
frequencies  for  a  source  at  larger  depths. 

Figures  7  and  8  contain  also  the  phase  spectrum.  It 
should  be  noted  that  troughs  in  the  amplitude  spectrum 
usually  coincide  with  irregular  variations  in  the  phase 
spectrum. 


i- 


hos . 
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4 .4  The  Influence  of  ndnor  variations  of  the  crustal 
model  “ 

To  get  an  estimate  on  how  minor'  changes  of  the  crustal 
model  might  influence  the  transfer  function,  variations  have 
been  introduced  into  the  standard  model  NOCR  (see  table  l). 

The  models  LVCl  and  LVC2  are  given  in  tables  2  and  3. 

Table  2:  Crustal  Model  LVCI 


Depth  interval 

(^0 _ 

0-3 

.3-6 

6  -11 

11  -  20 

20  -  45 

45-00 


P- 

'  km/sec 
3.10 
5.95 
5.50 
6.03 
6.97 
8.10 


S-Velocity 
(km/sec ) 

1.79 

3.44 

3.30 

3.48 

4.02 

4.67 


Density 


2.35 
2  65 
2.70 
2.60 
3.00 
3*35 


Table  3:  Crustal  Model  LVC2 


Depth  interval 
_ 

0  -  3 
3  -  6 
6-11 
11  -  20 
20  -  45 
45  '■  OO 


p- 

km/sec ) 

S-Veloclty 
(km/sec ) 

Density 

(^cm3) 

3.10 

1.79 

2.35 

5.95 

3.44 

2.65 

5.00 

2.90 

2.50 

6.03 

3.48 

2.60 

6.97 

4.02 

3.00 

6.10 

4,67 

3.35 

f  velocity 

channel  has  been  introduced 

The  velocity  in  the  LVC 

2  channel 

channel 

thus  increasing 

the  veloc- 

ity  contrast.  In  figure  9  a  comparison  between  the  transfer 
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functions  of  the  models  NOCR,  LVCl  and  LVC2  in  the  frequency 

range  0  to  0.3  cps  is  presented.  The  source  has  been  kept 

at  a  depth  of  8.0  km  for  the  three  models  so  that  for  LVCl 

and  LVC2  the  source  is  located  in  the  low  velocity  channel. 

0 

The  angle  of  incidence  is  Jp  =30  .  F’igure  9  merely  demon¬ 
strates  that  in  this  range  of  frequency  the  minor  variations 
so  introduced  in  the  model  are  virtually  of  no  effect  on  the 
transfer  function.  Closer  examination  shows  small  differ¬ 
ences  which,  however,  would  have  no  significance  in  signal 
analysis.  .In  the  frequency  range  greater  than  about  1  cps 
differences  in  the  transfer  function  begin  to  become  notice¬ 
able.  Prom  the  analysis  of  these  and  other  numerical  exam¬ 
ples  the  following  might  be  concluded  about  the  influence  of 
minor  variations  in  tne  crust: 

Minor  variations  in  the  crust  become  significantly 
apparent  in  the  transfer  function  only: 

1)  if  the  velocity  contrast  of  the  newly  Introduced 
layer  is  sufficiently  large. 

2)  if  l)  is  fulfilled,  only  that  part  of  the  spec¬ 
trum  will  be  affected  above  a  frequency  |  where 
^0  is  about  the  fundamental  frequency  of  "destruc¬ 
tive"  Interference  between  reflected  waves  in  the 
new  layer. 

5.0  Conclusions 

Analytical  expressions  for  the  transfer  functions 
for  dilatations 1  body  waves  radiating  into  the  mantle  from 
a.,  system  consisting  of  a  point  source  in  a  layered  crust  have 
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been  established  for*  large  distances.  Three  types  of  point 
sources  have  been  treated:  an  explosive  source,  a  single 
couple,  and  a  double  couple.  • 

Preliminary  numerical  calculations  of  the  transfer  func¬ 
tion  for  explosive  point  sources  have  been  presented.  They 
show  as  esp^^cially  important  features  the  strong  influence 
of  the  source  depth  on  the  shape  of  the  spectrum  and  the  re¬ 
latively  strong  rejection  of  low  frequencies  for  shallow 
explosions . 

Further  analysis  of  other  models,  especially  of  the 
other  two  point  source  types  must  show  to  what  extent  the 
present  conclusions  may  be  generalized. 
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